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Abstract 

In tliis paper, we study a macroscopic system of electrically interacting metallic beads organized 
as a sequence along an annulus. A random mechanical shaking mimics the thermal excitation. We 
exhibit non Fickian diffusion (Single File Diffusion) at large time. We measure the mobility of the 
particles, and compare it to theoretical expectations. We show that our system cannot be accurately 
described by theories assuming only hard sphere interactions. Its behavior is qualitatively described 
by a theory extended to more realistic potentials [Kollmann, PRL 90 180602, (2003)]. A correct 
quantitative agreement is shown, and we interpret the discrepancies by the violation of a key 
assumption of the theory, that of overdamped dynamics. We recast previous results on colloids 
with known interaction potentials, and compare them quantitatively to the theory. Focusing on 
the transition between ordinary and single file diffusion, we exhibit a dimensionless crossover time 
that is of order one both for colloids and our system, although the time and length scales differ by 
several orders of magnitude. 

PACS : 05. 40. -a Fluctuation phenomena, random processes, noise, and Brownian motion 66. 10. eg 
Mass diffusion, including self-diffusion, mutual diffusion, tracer diffusion, etc. 47.57.eb Diffusion 
and aggregation 
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I. INTRODUCTION 



When particles transport occurs in highly confined geometries, in such a way that they 
are distributed along a line, and contained in channels so narrow that they cannot cross 
each others, very strong correlations appear between the particles, and diffusive transport is 
anomalous. The particles are said to undergo Single File Diffusion (SFD). Such occurences 
include transport of water and ions through molecular-sized channels in biological mem- 
branes 1-3], molecular sieving effects in nanoporous zeolites j4-6|, diffusion of molecules 
along folded polymer chains j^l and diffusion of colloids in microfluidic devices js, 9| or in 
specially designed channels generated by scanning optical tweezers [iG 



Hi. 



SFD is a long standing problem in probability theory. The simplest modelisation con- 
sider the particles as hard spheres, interacting with a zero range singular potential 12N21|. 



and revelled the prominent effect of the confinement : The diffusion of a given particle is 
subdiffusive, with the mean-squared displacement (Ax^) that scales at large times t with a 
non-fickian exponent, 

(Ax2) = Ff/\ (1) 



F beeing called the mobility. For hard core interactions between particles, the mobility 
reads 

F« - ^,/^. (2) 



2 Do 

where p is the number of particles per unit length and Dq the single particle diffusivity. 
This non-fickian scaling at large time has already been observed for molecular transport 
{4-6], colloids ine one-dimensional channels 



in zeolites (4-l6|, col 
in circular channels 
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-Qjj, charged macroscopic beads 



23| and in numerical simulations of brownian particles interacting 



with a finite range pair potential 24j-l26| . When the focus is put on the mobility itself, 
some of those works show discrepancies with the hard sphere description ([2]). Gupta et al. 
measured the mobility of ethane in AIPO4 — 5, and found it one order of magnitude greater 



than expected 



J]. For colloids, the measurements show a dependence of the mobility on the 



26( 1 . in contrast with 



strength of the interparticle potential 

n 

Recently, KoUmann [27|] extended SFD theory to brownian particles interacting with 
a finite range pair potential. His work relies on three fundamental assumptions. First, 
the calculations are valid in the thermodynamic limit, for infinite systems. Then, at the 
very basis of the calculation, it is assumed that the dynamics of the brownian particles 



is overdamped. Lastly, the hydrodynamic interactions, which are the feedback effect on a 
given particle motion that comes from the flow of the surrounding fluid induced by that very 
motion, are neglected. In this paper, we want to address the relevance of those assumptions 
to actual instances of SFD. We discuss below finite size effects and overdamping, and we 
review previous works in which hydrodynamic interactions could be an issue |8|. 

Another issue is the dependence of SFD mobility on the temperature, brownian particles 
density and interaction strength. The evolution of the mobility with the density, for finite 
range interactions between the particles, have been checked for colloids by Lutz et al. 
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Hi, 



and by Lin et al. 9|, who found good agreement with KoUmann's calculation. It was 
pointed out in 9] that some colloids behave as a Tonks gas [2^, i.e. a gas of hard rods, 
for which Kollmann's result is equivalent to the classical hard sphere model, rescaling the 
density to take into account the size of the particles. Other authors have varied the strength 



of the interactions between the brownian particles in experiments 



and simulations 



26(1 . They observed the dependence of the mobility on the interaction strength but do not 



compare their results to Kollmann's calculation. We review their works, recasting their data 
to allow quantitative comparison with theoretical predictions. Then, we discuss our own 
experimental observations of SFD in a system of macroscopically charged particles in which 
the temperature, particle density and interaction strength may be varied independently. 
The interaction potential in our system is well described by a modified Bessel function Kq 
[see 291] and eqn. (fTOl) below], which is of the same functional shape as vortices interactions 
in superconductors. This study may thus be relevant to vortex transport in mesoscopic 
channels in superconductors 3Ch32|. 

The paper is organized in the following way. In section [Tll we review the theoretical pre- 
dictions for SFD with finite range potentials between the diffusing particles (§ 111 Ap . describe 
the relevant interaction potentials and recall their link with the isothermal compressibility 
(§ III Bl) . and compare the observations to the predictions of Kollmann in III Ci Then, in 
section UTTl we describe our experimental system (§ 1111 Al) . insisting on the random mechani- 
cal shaking applied to the system, that behaves as an effective thermodynamic temperature 
(§ nil B|) . In section HVt we set out our results, providing evidences of SFD for macroscopic 
charged beads (§ IIV AO . and studying the evolution of the mobility with the strength of in- 
terparticle potential, particle density and (effective) thermodynamic temperature (§ IIVBI) . 
We discuss potential finite size effects (§ lIVCp and put peculiar emphasis on the crossover 
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time between normal diffusion and SFD (§ lIVDp . A last section |V] is devoted to a summary 
of our results and our conclusions. 

II. SINGLE FILE DIFFUSION : A STATE OF THE ART. 
A. Single file diffusion for finite range interparticle potential 

Recently, KoUmann j27| have extended SFD to more realistic, finite range interactions 
between the diffusing particles. The subdiffusive behavior, with a mean-squared displace- 
ment (hereafter "m.s.d. ") scaling at large times t with a non Fickian exponent 1/2 as in 
([2]), is recovered. However, the mobility is now given by 



F, = -5(0,0)J^, (3) 
p V vr 

where we have introduced the static structure factor of the particles, in the limit of long 
wavelength or zero wavenumber q, S{0,0) = S{q — )■ 0,t = 0). The diffusivity -De// is the 
effective diffusivity of a brownian particle, taking into account its interactions with the other 
brownian particles, and differs from the single particle diffusivity -Do- Independently from 
peculiar effects linked to SFD geometry, interparticle interactions lead to a renormalization 



of the single particle diffusivity 



33] . When hydrodynamic interactions are neglected, the 



effective diffusivity -Dg// = Dq/S (0,0) 33|. With this expression, the mobility of brownian 



particles with finite range interactions reads 



p V vr 

The static structure factor may be measured by light diffusion for suspensions, or calcu- 
lated from the particles positions in simulations. It may also be expressed in terms of the 
isothermal compressibility kt as S'(0,0) = phsTuT, where /c^ is the Boltzmann constant 



and T the thermodynamic temperature 33|. This is particularly well suited to the analysis 
of our experiments, since we have rather small systems and a poor resolution on structure 
factor measurements. In contrast, the interaction potential is known with a good precision, 
and allows an easy determination of kt- 

Generically, the interaction potential may be expressed as f/(r) = Uof{r/\) where A is 
some characteristic length, Uq some energy scale and r the distance between the particles. 
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For interacting particles on a line with a mean interparticle distance (r) = 1/p, we may 



calculate the isothermal compressibility at zero temperature 



p 



(5) 



m>0 ^ m>0 

where the summation extends, in principle, to all particules. In practice, for the potentials 
that we consider, the sum converges fastly and keeping only three terms in the sum already 
gives four significant figures. The mobility thus reads 

n 1/2 



P 



2 DopkBTKT 



TT 



D^kBTX' 



m>0 



(6) 



This result is valid for a suspension of brownian particles interacting with a pair potential 
that is a function of interparticles distance only, when hydrodynamic interactions due to 
the coherent solvent motions are neglected. In principle, this calculation does not include 
entropic contributions since it is valid at zero temperature, but finite temperature corrections 



are small 



35|. 



For the sake of comparison with other works l8|, |26i, we recall the definition of the 
reduced mobility F* = pF/ ^/Dq, which for hard sphere potential is a constant equal to 
©• We also introduce the dimensionless potential F = f3Uof{l/ Xp) where (3 = 
l/^ksT). With those variables, the reduced mobility may be writen 

2 pX 



,m>0 ^ 



If 



Xp 



1/2' 



(7) 



It is clear from this expression that, although dimensionless, the variables (F,F*) are not 
scaling variables since there remains a dependence of the reduced mobility on the density p. 



B. Interaction potentials 

For the sake of convenience, we recall here the interaction potentials of concern in what 
follows : 

(a) the interaction between magnetic moments induced in a brownian suspension of para- 
magnetic particles by an applied magnetic field B The dimensionless potential energy 
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Tpara reads 
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with /io is the vacuum permeabihty, Xeff the magnetic suceptibihty of the particles, and 
where the characteristic length A may be, e.g. the particle diameter a. The susceptibility 
that corresponds to the experiments |8| is 2.2 x 10^^^ A m^T~^, the same value beeing used 



in the simulations I26l. 



(b) a screened electrostatic interaction described by a DLVO pair potential (see 33| . 
appendix A), with dimensionless potential energy Tdlvo 



2 Xb [ e'^/^^A) \ 2 g-r/A 
DLVO = ^e//— ^YT^7(2A)j 



where Ze// is the effective charge of the colloid, Xb the Bjerrum length, A the Debye screening 
length and a the particles diameter. The relevant values of those parameters, used in the 



simulations 



24| and taken from experiments [36|, are a = 2.8 /im, Z^ff = 5400, A^ 



0.72 nm and A = 550 nm. 

(c) the interaction energy for our experimental system of metallic beads contained in a 
condenser of thickness h under an applied voltage Vq. The set-up is described below in 
sec. IIII A[ The dimensionless pair potential Fbeads reads |29| , 

Tbcads = ^eoeoV.'hKo {£) , (10) 

where cq ~ 0.71 is a numerical constant, eo the vacuum permittivity, Kq the modified 
Bessel function of second order. In our experiments h = 1.5 mm and the characteristic 
length A ~ 0.32/i ^ 0.5 mm. The parameters eo and A depend on the dimensions of the 
beads, and the geometry of the condenser but are independent of the applied voltage Vq 
(see below sec. IIII Ap . They are kept fixed in all our experiments. The tension Vq rnay 
vary in the range [800, 1400] V, allowing us to easily monitor the interaction energy between 
the metallic beads. This interaction is similar to the screened electrical interaction between 
charged particles in suspensions ([9]), because it depends on r as exp{—r/X)/y/r for large r. 



C. Experiments and simulations on SFD with soft interactions 



In this section, we discuss previous works on colloids in SFD geometry, and compare 

27] . We have collected the 



quantitatively their results to the predictions of KoUmann 



6 



existing data and plotted t 



le measured values of the reduced mobility F* [see ([7])], either for 



24j . In all cases, the interaction 



varying potential strength [8|, |26[ or varying packing fraction 
potentials are explicitely known (see § 111 Bp and equation ([7]) allows the calculation of the 
reduced mobility predicted by Kollmann without any free parameter. We add to the existing 
mobility data the corresponding theoretical predictions, thus providing a quantitative test 
of the theory. We recall in Table [T] the relevant experimental or numerical characteristics of 
the systems that are discussed below and in section IIVI 



System 


Reference 


U{r) 


N 


Dynamics 


H.I. 


Theory 




finite range 


oo 


overdamped 


no 


Colloids (num) 


m 


m 


900 


over damped 


no 


Colloids (exp) 


[8] 


m 


89 


overdamped 


yes 


Colloids (num) 


[26] 


m 


900 


overdamped 


no 


Colloids (exp) 


[10, 


m 


49 


overdamped 


no 


Colloids (num) 


[25] 


m 


G [50, 100] 


underdamped 


no 


Charged beads (exp) 


this work 


m 


G [12,37] 


underdamped 


no 



TABLE I: Characteristics of experimental (exp) and numerical (num) systems undergoing SFD. 



Colloid dynamics is presumably overdamped in experiments 
Interactions" . 



33| . H.I. stands for "Hydrodynamic 



1. Colloids with paramagnetic interactions 

Wei et al. realized single file geometries by putting colloids in circular channels etched by 
photolithography. The colloidal particles are polystyrene beads doped with Fe203 clusters, 
so that when submitted to an external magnetic field B they interact as magnetic dipoles 
whose moments depend on i?, with an interaction potential given by (|8]). It is thus possible 
to study the behavior of the system with varying interaction strength. The same system is 



simulated by Herrera- Velarde and Castaneda-Priego [26 1. 

Single file diffusion, with a non-Fickian scaling of the m.s.d. as t^^"^ is observed in both 
experiments and numerical simulations, allowing the determination of the mobility. A de- 
pendence of the reduced mobility as a function of the reduced interaction Fpara is observed 



for both data sets (see Fig. 4 of 8| and Fig. 3-b of 26|]), summarized here in our Fig. [T] 

0.4 r 




FIG. 1: (color online) Plot of the reduced mobility F| as a function of the reduced interaction F 



for colloids with magnetic dipole interactions. • : Numerical simulations 
data jsl, 



26l |. ■ : experimental 



solid line Kollmann's formula, Eqn. ([7]) with ([8|) (no free parameter in this calculation) 



The numerical simulations fulfill all assumptions of the theoretical calculations (see Ta- 
ble HI), and their results are indeed in excellent agreement with the theory. The small 
discrepancy between experiments and theory may be due to the hydrodynamic interactions 
between the particles in the file. The experiments are done at constant density and constant 
temperature, hence the fact that the agreement is less good at small magnetic potential 
(rpara < 4, typically) may be explained by the greater relative importance of hydrodynamic 
interactions compared to those between magnetic dipoles. Hydrodynamic interactions are 



basically repulsive, so that neglecting them lead to overestimating the compressibility o 
brownian particles. This is consistent with the fact that the SFD mobility measured in 



the 
is 



less than what is expectated without hydrodynamic interactions. In the experiments of Lutz 



et al. 



11| , great care have been taken to confine the beads in optically designed channels, in 



order to minimize hydrodynamic interactions. The quantitative agreement with Kollmann's 
results seems better in this configuration. 



2. Screened electrostatic interactions 



In another set of numerical simulations 2^ , Herrera- Velarde and Castaneda-Priego study 
colloidal particles with screened electrostatic interactions, described by the DLVO pair po- 
tential ([9]). They vary the packing fraction ip of the brownian particles, which is linked to 
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the density phy p = (f/a, where a is the particles diameter. They observe SFD behavior, 
with a non-fickian scahng for the m.s.d. (see their Fig. 2-a), and measure the SFD reduced 
mobihty ([7]) as a function of the colloid packing fraction (p. 
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FIG. 2: (color online) Plot of the reduced mobility Fg as a function of the brownian particles 
packing fraction for colloids with screened electrostatic interactions (DLVO). • : Numerical 



simulations 
lation) . 



24( 1 ■ solid line Kollmann's formula, Eqn. d?]) with ([9]) (no free parameter in this calcu- 



In Fig. O we compare the simulations data (see Fig. 2-b of |26|) to the theoretical pre- 
diction of KoUmann. The agreement is very good at high packing fraction {ip > 0.4), and 
very poor at lower ones. A possible issue explaining the discrepancy at low packing fraction, 
is that the time taken to fully developed non-fickian scaling for the m.s.d. may be 

large compared to the duration of a simulation. When particles are organized sequentially 
on a line, their diffusion is normal (Fickian) at small time, and they must cover diffusively 
a lengthscale comparable to their mean distance in order to feel the presence of other beads 
and undergo SFD. SFD cannot be observed before a typical time of order a'^/D^ff where a 
is the interparticle mean distance and -De// the diffusivity of interacting particles (see below 
section HV Dp . In simulations as in experiments, it means that the measurement duration has 
to be greater at smaller packing fraction. Indeed, as seen from Fig. 2-a of 2J|, at a packing 
fraction of 0.288 the t^^^ scaling is not very good, and the reduced mobility measurement 
may not be very good either. This effect should be worse at lower packing fraction. 

Involving the smallness of the packing fraction is consistent with other results. The results 
reported in Fig. [1] are obtained with paramagnetic colloidal particles, at a packing fractions 
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of 0.33. For our own system, the packing fraction is always greater than 0.38. As shown by 
Fig. [6] or Fig. [71 the agreement with Kollmann's result is much better than for simulations 
of low packing fraction colloids (Fig. 



3. Discussion 



We have done in this section a quantitative comparison between experimental and nu- 



merical studies of systems involving SFD, anc 



a theory that takes into account finite range 



interactions between the diffusing particles j27|. The qualitative behaviors, predicted by the 
theoretical analysis, are indeed observed : At large time, non-fickian diffusion occurs with 
the peculiar SFD subdiffusive behavior (Ax^) = Ft^^"^ for the mean squared position fluc- 
tuations. The mobility F decreases when the brownian particles compressibility decreases, 
either because of an increasing interaction strength (§ 111 C II) or of an increasing packing 
fraction (§ III C 21) . This shows that the assumptions at the basis of the theoretical approach 
are rather consistent with actual realizations of SFD (see Table [I]). 

Quantitatively, the test is stringent because once the potential is known, then the re- 
duced mobility is determined without any free parameter. For colloidal particles interacting 
as magnetic dipoles js], at a packing fraction of 0.33, the agreement of theory with 
simulations is excellent, and very good with the experiments (see Fig. [T]). For screened elec- 



trostatic (DLVO) interactions, the simulations are done with varying packing fraction 24j . 
and the agreement is good at high packing fraction only (see Fig. |2]). 



In the next sections, we do this quantitative comparison for a versatile system (see IIII Ap 



of macroscopic charged beads in which the interaction strength, temperature, density and 
system size may be independently varied. The interaction potential is known from previous 
studies 29j, and the comparison again involves no free parameter. The time resolution and 
experiment duration is sufficient to describe clearly both the short time regime, intermediate 
between ballistic transport and normal (Fickian) diffusion, and the long time regime for 
which we exhibit SFD behavior. We discuss the crossover time between those regimes, 
and show that, in order of magnitude, it is the time taken by a tagged particle to explore 
diffusively the potential well created by its two neighbors. When the SFD is observed, 
we measure the mobility and discuss its variations with the potential strength, the density 
and the temperature. We pay particular attention to potential finite size effects and show 
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that, even for quite small systems, they do not intervene because of the periodic boundary 
conditions. This is a solid evidence that calculations done in the thermodynamic limit, hence 
infinite system size, are actually relevant. 



III. A SYSTEM OF MACROSCOPICALLY CHARGED PARTICLES 
A. Experimental set-up 

Our experimental set-up consists in millimetric stainless-steel beads (diameter 0.8 mm, 
mass m = 2.15 mg) located on a silicon wafer which is the bottom electrode of a plane 
horizontal condenser. A metallic frame intercalated between the electrodes, in contact with 
the bottom one, confines the beads in a circular channel of radius R in such a way that they 
cannot cross each others. The bead motions may be decomposed in radial and orthoradial 



coordinates, wich have been shown to be mutually statistically independent [22(]. Since the 
annulus radius of curvature is much greater than the beads mean distance, we have single file 
diffusion with periodic boundary conditions. We have done experiments with two annulus, 
one radius R = 4 mm and width 2 mm, the second of radius R = 9 mm and width 2 mm. 
We report on experiments done with N = 12, 14 or 16 beads in the small annulus (densities 
p = 477, 557 and 637 particles per meter), and with = 27, 32 or 37 beads in the large 
one (p = 477, 566 and 654 particles per meter). 

When a voltage Vq G [0.5, 1.3] kV is applied between the two electrodes, the beads 



become charged and repel each other with a known interaction [29| [see above eqn. ( ITOj) ]. 
In order to introduce an effective thermodynamic temperature, the condenser is fixed on 
loudspeakers excited with a white noise voltage. We have verified that the mechanical 
shaking, transmitted to the beads via the friction with the bottom electrode, does indeed 
behave as a thermal noise. This point is discussed in the following section IIII Bl and in a 
preceding work j37|. 

To avoid any gravity effects on the beads, we have to ensure precisely the horizontality of 
the annulus. The experimental cell is supported by three columns, the extremities of which 
define a plane, the vibrations that mimic thermal excitation are exerted in that same plane, 
and we ensure a priori the horizontality of this plane. However, no inclination measurement 
is more sensitive than the system motion itself : the mean particle distance should be uniform 
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along the annulus for a satisfactory horizontality, which is checked before each experimental 
run. The picture displayed in Fig. [3] shows such a situation. 

The top electrode is a glass plate covered with an ITO (Indium-Tin-Oxyde) layer which 
is sufficiently fine (thickness 0.1 fim) to ensure optical transparency. Images of the particles 
may thus be recorded in real time during the experiments. A typical snapshot with 37 beads 
is provided in Fig. |3l The interval between two snapshots ranges between 10 and 30 ms and 
series of 10 000 pictures were recorded. The effective thermal bath is characterized by a 
relaxation time about 100 ms 22|, hence individual trajectories of the beads are determined 
over long times. The diffusion of the beads in an annulus of radius R is measured through 
the time evolution of their orthoradial m.s.d. R^{A6'^{t)) given by 

R^Ae^it)) = R\ [9{t + to) - e{u) - {e{t + to) - eiu))] ') (ii) 

where 9 is the orthoradial cumulated angle in radians and to an arbitrary initial time. The 
bracket (■) denote ensemble averaging. This averaging is done on every beads, since they 
all play an equivalent role. Moreover, the phenomenon is assumed to be stationary, so that 
{A6^{t)) do not depend on to. For a given time t, it makes thus sense to average on the 
initial time to. Let n be the overall number of data recorded in one experimental run, 6t the 
sampling time, and Ut = t/6t. Then the averaging on the initial time to reads 

^ n-nt + l n-rit + l /' 

1=0 \ 1=0 / 

where the index i is such that to = iSt. This way of averaging greatly improves the statistics 
when Ut is smaller than n. In what follows, this double averaging will be denoted {A6^{t)) 
for simplicity. 



B. The effective temperature 



It has been carefully checked in previous works 22|, |37| that the mechanical agitation 
is, indeed, equivalent to a thermodynamical temperature. In 37|, the random position of a 
single bead rolling on a wafer inclined by a very small angle was recorded, and its probability 
distribution compared to the Boltzmann distribution for several shaking amplitudes A. An 
excellent agreement was obtained, allowing the definition of an effective temperature T 
which depends on the shaking amplitude. Then, an in situ thermometer was built, made 
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FIG. 3: Picture of the circular channel with = 37 beads inside, recorded during an experimental 
run (T = 10.2 10^^ K, Vq = 1000 V). Inner radius 8 mm, outer radius 10 mm, beads diameter 
0.8 mm. In the right-up corner, an identical bead confined in a disk of radius 5 mm is the embedded 
thermometer. 

from a single bead trapped in a circular frame located near the main one and submitted to 
the same voltage. Such a thermometer is shown in Fig. [31 The bead is confined in a 2D 
parabolic potential. As it should, its radial mean square displacement is a linear function of 
the effective temperature T in the long time limit, hence this set-up actually behaves as a 
thermometer. The temperature calibration is shown in Fig. |4](a). The voltage Vq = 1000 V 
was kept constant during all measurements. The error bars (roughly ±5% of the mean value) 
reflect the statistical dispersion of the measured value for several experiments undertaken 
at the same shaking amplitude A. 

The mobility of a tagged particle in a file depends on the single particle diffusivity Dq, 
both for hard spheres ([2]) or finite range interactions. In this latter case, the effective 



diffusivity -De// is a function of Dq, the temperature and the compressibility (jlj). In 22[, the 
diffusion of a single bead in a circular channel was studied. It was shown that the radial and 
orthoradial motions are completely decoupled. The orthoradial motion is well described by 
a free diffusion with diffusivity Dq, whereas the radial motion is that of a brownian particle 
in a ID parabolic potential. The diffusivity Dq has been measured as a function of the 
mechanical shaking amplitude A, at constant applied voltage. As shown by Fig. H] (b), it is 
a monotonously increasing function of A. 

The mechanical shaking exerted on the beads acts as an effective thermal bath, and it is 



thus possible to describe the beads motion with the help of a Langevin equation (22|. The 
complicated coupling, due to solid friction, between the horizontal vibrations of the silicon 
wafer and the bead motions behaves as the solvent particules that exert brownian random 
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FIG. 4: (a) Effective thermodynamic temperature T (in 10^^ K) as a function of mechanical 
agitation A (a.u.). (b) Single particle diffusivity Dq (10~® m^ s~^) for orthoradial motions in 
an annulus, as a function of mechanical shaking amplitude A (a.u.). In each graph, error bars 
are determined statistically from several measurements, and represent typically ±5% of the mean 
value. 

forces on a colloid. However, the shaking is physically different from the erractic motion of 
solvent molecules around colloidal particles, and hydrodynamic interactions are avoided in 
our system. 



IV. SINGLE FILE DIFFUSION OF MACROSCOPIC CHARGED PARTICLES 



A. Single file diffusion 

In Fig. [5l we plot the orthoradial m.s.d. R'^{A6^) as a function of time, for p = 654 m~^, 
T = 10^^ K, and varying applied voltage Vq, hence varying interaction strength [see ffTOj) ]. 
The behavior of a bead in a single file geometry depends on the observation time. At very 
small time, the bead motion is ballistic, with a m.s.d. oc t^. At intermediate times, normal 
(Fickian) diffusion with m.s.d. oc t takes place, because some time is needed for the statistics 
of the particle to be influenced by the single file geometry. As shown by the inset, in this 
regime the m.s.d. is somewhat faster than linear, intermediate between ballistic transport 
and ordinary diffusion. This behavior is independent of the applied voltage. 

For brownian particles confined on a line, SFD may occur with non-fickian scaling 
m.s.d. oc t^/^ for asymptotically long times. This is what we observe at times greater 
than typically 200 ms, as shown by Fig. [5l for all voltage values. The inset gives evidences 
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FIG. 5: (Color online) Plot of the orthoradial m.s.d. i?^(A^^) (in mm^) as a function of time t 
(in ms), in a log- log scale, for p = 654 m"^ (37 beads, R = 9 mm), and T = 10^^ k. The applied 
voltage Vo vary from 900 V to 1400 V in steps of 100 V in order blue-black-red-green-pink-orange. 
For clarity, each data set have been shifted by Log(2) except for Vq = 1400 V. The solid lines are 
the best fits with ([T|), where the mobility F is the adjustable parameter. The inset display the 
unshifted data, showing that the diffusion at small time is independent on the voltage. The solid 
line is of slope 1. The dimensionless potential energy is such that 16 < Fbeads < 40. 

that the mobility actually depends on the applied voltage, hence on the interparticle inter- 
actions. This dependence is studied in detail in section IIVBI The main plot shows that the 
quality of the fit is good. It evidences a crossover between a normal diffusion regime, and 
the appearance of SFD behavior. We may thus define a crossover time Tc, which is shown 
to increase with decreasing applied voltage, hence decreasing interparticles interaction [see 
( fTOj) ] or increasing compressibility [see (|5])]. The same behavior was already observed in 
colloids (see jsj. Fig. 2-B). A more quantitative analysis of the crossover time is the subject 
of section IIVDI below. 



B. Mobility measurements 

The t^^^ scaling of the m.s.d. , typical of SFD, is clearly observed in Fig|5l allowing 
mobility measurements using ([T]) with the mobility as the fitting parameter. We vary the 
applied voltage and the mechanical agitation, and to increase the statistics, we record several 
picture sets at the same voltage and shaking amplitude. We display in Fig. IHlthe mobility as 
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a function of the applied voltage, for p = 637 [Fig. [6]-(a)] and p = 477 [Fig. |6]-(b)], 
at fixed shaking amplitude. As shown by the error bars in Fig. HI there is some uncertainty 
on the temperature for a given mechanical agitation level. It means that, when we want 
to compare our mobility measurements to the theoretical prediction (jl]), the error done 
when evaluating the temperature T and diffusivity Dq, knowing the shaking amplitude A, 
prevails on any other one. We thus plot the prediction from for the lowest and highest 
temperatures consistent with amplitude A. 

The data clearly exhibit qualitative features that are in excellent agreement with the 
theoretical predictions, at both densities. The mobility depends on the strength of the 
interparticle interaction, showing that SFD cannot be described by the mobility ([2]), valid 
for hard sphere interactions. The mobility decreases when the voltage increases. This is 
fully consistent with because increasing Vq means increasing the interaction energy and 
decreasing the compressibility, as shown by (fTO!) and ([5]), hence decreasing the mobility as 
stated by (j4j). Quantitatively, the agreement is correct. An important point, that will be 
developed in the next section, is that the system is described by its density rather than by 
its size, because the mobility values taken with small and large annulus are consistent at a 
given density. 

In Fig. [TJ we display the mobility as a function of temperature, for a constant density 
of 561 m~^, and two voltages. Since the diffusivity Dq{T) is an increasing function of the 
temperature, the mobility should increase with the temperature, and it should be higher 
at lower applied voltage since the compressibility decreases with the voltage. At a given 
temperature, the mobility measured for Vq = 1300 V is indeed always smaller than for 
Vq = 1000 V, in good qualitative agreement with the theoretical picture. However, the 
quantitative agreement is not very good. 

In order to interpret those discrepancies, two possible issues may be considered. The first 
one is the smallness of the beads number, whereas calculations are done in the thermody- 
namic limit. The larger beads numbers in our experiments (27, 32 and 37) are nevertheless 
quite similar with those in previous experiments (see Table [T]). 



Another issue is suggested by Nelissen et al. [25| to interpret the discrepancies with SFD 
behavior that they observe in their simulations. Their algorithm use the Langevin equation 
with DLVO interaction between particles without assuming overdamped dynamics, taking 
into account the inertial effects (acceleration). Overdamped dynamics is at the basis of 
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FIG. 6: (a) Mobility Fs, in mm^s~^/^ as a function of the applied potential Vq in V. Blue dots, 
N = 37, R = 9 mm; red dots, = 16, i? = 4 mm. The solid lines are the theoretical expectations 
for temperatures T = 10.6 10^^ K and T = 11.7 10^^ K (no free parameter). Dimensionless 
potential energy : 11 < Fbeads < 55. (b) Blue dots, = 27, i? = 9 mm; red dots, N = 12, 
i? = 4 mm. The solid lines are the theoretical expectations for temperatures T = 10.7 10^^ K and 
T = 10.2 10^^ K (no free parameter). Dimensionless potential energy : 6 < Fbeads < 10. 



26|. It is a safe assumption 



both Kollmann's calculation, and the algorithm used in 
in colloids 33[, but our system is clearly not in this regime, as shown in 22|. When the 
radial m.s.d. of a bead in a parabolic potential is compared to the predictions of a Langevin 



equation (see the Fig. 4 and Fig. 5 of |22|), it is clear that the bead undergoes oscillations, 
that are not overdamped since a pseudoperiod should easily been defined. We may get an 
intuitive idea of the corrections due to inertial effects as follows. If the particle motion is 
not overdamped, it declines more slowly than in an overdamped regime. It may explain why 
we always observe a greater mobility than the prediction obtained assuming overdamping. 
Another satisfactory issue is that the agreement is better either at low temperature (the 
lowest temperature in Fig. [7] is in very good agreement with the theory) or at large applied 
voltage (hence large interaction), when thermal motion is comparatively less important. 



C. Finite size effects 



The theorerical description of SFD 27|] takes place in the thermodynamic limit, for infinite 
particle number A^, system size L, and finite mean density p = N/L. In either experiments 
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FIG. 7: Mobility Fs, in mm^s'"^^^ as a function of the temperature (±5% error bar), for a density 
of p « 561m^^. Blue error bars : Vq = 1300 V and 11 < Fbcads < 19, red error bars : Vq = 1000 V 
and 6 < Fbcads < 11- The solid lines follow the same color code, and are the theoretical expectations 
for the mean temperature. 



or simulations, the particle number is obviously finite, and boundary effects are avoided 
using periodic boundary conditions, putting the diffusing particles in an annulus. However, 
a finite size effect that cannot be avoided is the fact that an acoustic mode propagating 
along a ring, comes back to its point of departure after one turn rather than going away 
at infinity. Such modes, obviously missing in infinite systems, may have an effect on the 
correlations that are responsible for the SFD. It is thus of interest to look at the behavior 
of systems with the same density, but different sizes. 

In Fig. ini we display mobility measurements done at two densities, both with large and 
small annulus. The smaller system contains few beads (12 and 16), but the data are nev- 
ertheless in good agreement with those for larger bead numbers and the same density. The 
density is thus the relevant parameter, and finite size effects are not exhibited by the mea- 
surements. In order to test potential finite size effects on more sensitive data, we display in 
Fig. [8] the rough recording of orthoradial m.s.d. as a function of time. All data are obtained 
at the same temperature, as measured by the embedded thermometer, same applied voltage 
and two different densities : 477 beads per meter with = 12 in the small annulus and 

= 27 in the large one, 637 beads per meter with A^ = 16 in the small annulus and A^ = 37 
in the large one. 
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FIG. 8: (Color online) Plot of the orthoradial m.s.d. B?{/S.6'^) (in mm^) as a function of time t 
(in ms), in a log-log scale, for densities 477 m~^ (blue dots, small annulus, N = 12; black dots, 
large annulus, N = 27) and 637 m"^ (red dots, small annulus, = 16; green dots, large annulus, 
N = 37). The solid black lines are of slopes 2, 1 and 1/2. The temperature is the same in all 
experiments (T = 11.2 10^^ K) and the applied voltage is 1000 V. The crossover time decreases 
when the density increases. 

The time evolution of the m.s.d. is the same for a given density, whatever the particles 
number. At large time, the m.s.d. evidences SFD behavior in each case. As may be ex- 
pected, the systems with smaller density exhibit the larger m.s.d. . All four data sets are 
necessarily recorded separately, because they involve a complete dismantling and reassembly 
of the experimental cell. They thus show a very satisfactory reproducibility, and provide a 
convincing evidence that the system is correctly described by its density, without finite size 
effects. 

As was already the case for the data reported in figJSl the small time evolution of the 
m.s.d. is faster than linear, between ballistic transport and ordinary diffusion. At larger 
times, non-fickian SFD with t^/^ behavior is recovered. It is clear that the crossover time 
between normal diffusion and SFD increases with decreasing density. This is satisfactory, 
because a smaller density means a larger distance between the particles, hence a longer time 
to feel their neighbors and undergo SFD. 
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D. Crossover time 



In this section, we focus on the crossover time Tc at which the system evolves from 
ordinary diffusion toward SFD. Since this evolution is continuous, the crossover time is 
roughly defined, in terms of order of magnitude only. This crossover time is the time needed 
by a given particle to explore a lenghscale sufficiently large to "feel" that it is in a sequence, 
and that its diffusion is perturbed by the neighboring particles. The relevant diffusivity is 
not Dq, that of a single particle in a thermal bath, but rather the renormalized diffusivity 



Z^e/z which takes into account the interactions between brownian particles [33|. Neglecting 
hydrodynamic interactions between particles, we get -De// = -Do/5'(0,0) = Dq/ [pksT Ht)- 
We thus expect the crossover time to be such that 

Using the relationship between the interaction potential and the compressibility ([5]), we may 
rewrite this expression as 

We can thus give another interpretation of the occurence of the effective diffusivity in (fT3|l . 
The motion of a tagged particle may be described as a diffusion with diffusivity Dq, in the 
potential U (r) due to the interactions exerted by its two neighbors (neglecting the obvious 
fact that the walls of this potential well are themselves fluctuating). If the particle is in 
a thermal bath at temperature T, at large times its typical m.s.d. in this potential well is 
kBT/U". The crossover time Tc is such that the particle "knows" that it is in a sequence, 
hence that it has felt the walls of the potential well due to its neighbors. This is exactly 
what is expressed by ( fT4l) . 

The compressibility may be calculated from ([5]) and the relevant expressions for the 
interaction potentials (E]), ([H]) and ( |TOl) . Several papers provide evidences of a transition 
between normal diffusion and SFD, in such a way that a crossover time may be estimated 



in each case 



0, y, Q, 



24 



26| . Moreover, we deduce the crossover time for macroscopically 



charged beads from Fig. [51 Those estimates are summarized in Table HTl We point out that 
the dimensionless quantity TcP^-De// is in each case roughly of order one, although in their 
respective physical units, the numerical values for all quantities Tc, p and Dejf differ by 
several orders of magnitudes in colloids and in our experiments. 
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Reference 


Tc 




^0 






[8], Fig. 2(B) 


10 s 


11 fiUl 


0.036 ^m^s"^ 


3 10-3 


0.32 


m, Fig. 2(a) 




5.7 pm 




2 10^2 


1.16 


[26], Fig. 3(a) 




6.1 




4 10-3 


0.39 


Fig. El 


0.2 s 


1.6 mm 


0.02 mm^s"^ 


1.6 10-3 


0.44 



TABLE II: Crossover time for several SFD observations. For col 
in ll| (Fig. 2) and ^ (Fig. 2) are quite comparable to those of 



£ids, the observed crossover times 
, for roughly the same densities, 



temperature and particle sizes. 



In Fig. ini we display the evolution of the crossover time with the applied voltage. To this 
aim, we plot i?^(A6'^)/t as a function of the time t. At very small time, the transport is bal- 
listic, the m.s.d. scales as and this function increases as t. At large time, SFD occurs with 
a m.s.d. scaling as t^/^, so that the function decreases as When ordinary diffusion takes 

place, the function should exhibit a plateau. As shown by the inset of Fig. |9l this plateau is 
quite small. In fact, the function R'^{A6'^)/t rather display a maximum, and the position of 
this maximum gives a rough estimate of the crossover time. The crossover time decreases 
with the applied voltage, which is in qualitative agreement with (fl^ . Quantitatively, the 
observed variation is weaker than what may be expected. 




''0.080 
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1000 1100 1200 1300 1400 
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FIG. 9: (Color online) Crossover time tq (in s.) as a function of the applied voltage (in V.) In 
the inset we plot the m.s.d. divised by the time as a function of the time. The crossover time is 
estimated as the position of the maximum of this curve. 



21 



V. CONCLUSION 



In this paper we discuss quantitatively experimental evidences of SFD in the transport 
of macroscopic charged beads, interacting with a screened electrostatic potential. To ensure 
periodic boundary conditions, the beads are aligned in an annulus narrow enough to avoid 
any crossing. For a time resolution of 10 ms and experimental runs duration about 1 h, the 
system timescales are such that we have access to both short time dynamics, intermediate 
between ballistic transport and usual diffusion, and long time dynamics that evidences single 
file diffusion. 

In this regime, the orthoradial mean squared displacement is such that i?^(A^^) = Ft^/'^. 



For finite range interactions, it has been shown recent 
F depends on the brownian particles compressibility 



y by KoUmann that the mobility 



271 ]. When the interaction potential. 



hence the compressibility, is explicitely given, the mobility calculation does not involve any 
free parameter. All qualitative features displayed by our mobility measurements are in 
excellent agreement with the theory : The SFD mobility decreases whith the interaction 
potential (Fig. [H]), decreases with the brownian particles packing fraction, and increases 
with the temperature (Fig. [7]). Quantitatively, we observe a SFD mobility that is somewhat 
greater than the theoretical predictions. We interpret this discrepancy using an argument 
taken in j25|. Those authors point out that a key assumption of the theoretical calculation 
is that the particle dynamics is overdamped, which is not the case in our system. 

The theoretical calculations are done in the thermodynamic limit. All systems actually 
exhibiting SFD involve a finite number of particles, and periodic boundary conditions. In 
experiments, the particles number may be quite small. Using two different experimental cells, 
we show that the data recorded in small systems (12 to 16 particles) are consistent with those 
for larger systems (27 to 37 beads), and that experiments done at the same temperature are 
quite reproducible (Fig|8]). Since all systems in which SFD may be a practical issue {e.g., 
ion transport in biological membranes jl-3|) presumably involve few diffusing particles, this 
constitutes a satisfying evidence of the relevance of Kollmann's calculation. 

In our experiments, we clearly observe both the short time behavior of the system, where 
ballistic transport give way to ordinary (Fickian) diffusion, and the long time behavior where 
correlations led to SFD. We are thus able to define a crossover time Tc between those two 
regimes. This may be done, perhaps less convincingly, in colloids where the short time 
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regime is quite hard to exhibit. The relevant physical parameters differ by several orders of 
magnitude in our systems and in colloids. Nevertheless, the dimensionless quantity TcP^Dejf 
is of order one in each system, where the effective diffusivity -De// takes into account the 
interparticles interactions (see § llVDp . This is in good agreement with the intuitive picture 
of SFD, which takes place when a tagged particle have had the time to explore diffusively 
its environment, and feel the presence of the neighboring particles, placed at an average 
distance 1/p. 



We review previous numerical 2J, |26| and experimental 



data on colloids that evidence 



SFD behavior. We recast those data and provide a quantitative comparison with the theory 
for finite range interactions. All qualitative features predicted by the theory are observed 
(Fig. [Hand Fig. H]). Quantitatively, the agreement is excellent for not too small packing 
fraction, as shown by Fig. [1] The numerical algorithms fulfill all theoretical assumptions 
(see Table [T]), and hydrodynamic interactions may explain the small discrepancies between 
experiments js] and theory. At smaller packing fraction, the numerical data do not agree 
with theory (Fig. [2]), which is presumably due to an insufficient simulation duration. 



The calculations of Kollmann 



27| rely on two fundamental assumptions : overdamped 



dynamics and no hydrodynamic interactions. Those latter are either avoided or negligi- 
ble in the works reported here. However, in some practical situations (transport through 
molecular-sized channels in biological membranes, molecular sieving effects in nanoporous 
zeolites, diffusion of colloids in microfiuidic devices), hydrodynamic interactions may play a 
significant part in the dynamics. In the diffusion of molecules along folded polymer chains, 
it is rather the overdamping assumption that may be incorrect. We have shown that it 
should give observable corrections of the SFD mobility. Theoretical studies should thus be 
pursued in those directions, but as it states, the theory of Kollmann that takes into account 
the softness on the interparticle interactions constitutes a very significant improvement on 
previous approaches assuming hard sphere interactions only. 
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